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ABSTRACT 

Tidal disruption of dark matter halos around proto-globular clusters in a halo of a small galaxy is studied 
in the context of the hierarchical clustering scenario by using semi-cosmological N-body/SPH simulations 
assuming the standard cold dark matter model (fio = 1). Our analysis on formation and evolution of the galaxy 
and its substructures archives until z = 2.0. In such a high-redshift universe, the Einstein-de Sitter universe is 
still a good approximation for a recently favored A-dominated universe, and then our results does not depend on 
the choice of cosmology. In order to resolve small gravitationally-bound clumps around galaxies and consider 
radiative cooling below T = IQ^K, we adopt a fine mass resolution (;7Jsph = 1.12 x IO-'Mq). Because of the 
cooling, each clump immediately forms a 'core-halo' structure which consists of a baryonic core and a dark 
matter halo. The tidal force from the host galaxy mainly strips the dark matter halo from clumps and, as a 
result, theses clumps get dominated by baryons. Once a clump is captured by the host halo, its mass drastically 
decreases each pericenter passage. At z = 2, more than half of the clumps become baryon dominated systems 
(baryon mass/total mass > 0.5). Our results support the tidal evolution scenario of the formation of globular 
clusters and baryon dominated dwarf galaxies in the context of the cold dark matter universe. 
Subject headings: galaxies: formation — galaxies: halos — method: numerical — dark matter 



1. INTRODUCTION 

Structure formation in the cold dark matter (CDM) dom- 
inated universe is characterized by hierarchical clustering, 
where small objects collapse early, and merge into larger sys- 
tems, such as galaxies and clusters of galaxies. Therefore, the 
formation process of small-scale structure is crucial to under- 
standing the formation of larger objects. The tidal effect from 
larger objects could affect star formation history in each sub- 
halo (gravitationally-bound systems in galaxies) and the final 
structure of the host objects. 

A^-body simulation is a suitable method to follow the co- 
evolution of substructure and host halos. iKravtsov et alJ 
(|^04) studied the evolution of dark matter subhalos within 
galactic CDM halos using A^-body simulations without a gas 
component. They found that the tidal force is essential for 
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the evolution of substructure on galactic scale. Baryons, how- 
ever in general, have more concentrated distribution than dark 
matter owing to their dissipative nature, and hence evolution 
of these subhalos should be strongly influenced by baryonic 
processes. 

Quantitative behavior of tidal interactions between a 
three-component (stars, gas, and dark matter) substructure 
(hereafter clump) and a host halo can be directly stud- 
ied by numerical simulations. However, typical mass res- 
olutions in previous A^-body/SPH (Smoothed Particle Hy- 
drodynamics) simulations of galaxy formation is 10^ - 
IO^Mq in mass, ^ 1 kpc in length, and the num- 
ber of particl es in a halo i s typically 10^ - 10^ per 
galax y, (e.g., Katz & Gunnlig QUlNavarr o & BenzllQQUlKaJz 
1992| iSteinmetz. Muller 1994 ril995t iNavarro & Steinme^z 
1997 i Weil. Eke & EfstathiQuI Il998t (Sommer-Larsen et al 



1999 



21M 



Steinmetz & Nava rrolll999t lEke. E fstathiou & Wriah 
Koda et al.ll2000aibt lAhadi et al]]2003albt iMeza et al 
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20031 ISommer-Larsen et alj 120031 lOkamoto et a ll l2003t 
^overnato et al. 2004; Robertson et al. 2004; Okamo to et alJ 
2005i) . This mass range is almost the same or larger than 
the masses of the first collapsed objects (~ IQ^Mq) (e.g., 
iTegmark et al. 1997; Yoshida et al. 2003). Therefore, these 
previous simulations could not address the hierarchical clus- 
tering processes from the first collapsed objects in the CDM 
universe. It is also not clear how globular cluster-size objects 
are formed and evolved in forming galaxies. 

In this paper, we study the evolution of clumps in the halo 
of a host dwarf galaxy with high resolution A^-body/SPH sim- 
ulations. Our cosmologically-motivated simulations are im- 
proved from previous simulations in terms of (1) particles 
masses, which are IO^Mq for SPH particles and ^ IO'^Mq 
for dark matter particles (cf. ^ lO^M© for SPH particles in 
^badi et al. 2003a), (2) inclusion of the low temperature ra- 
diative cooling below IQ'^K, which allows formation of cold 
clouds. Consequently, we can resolve structures larger than 
Mjeans > 2A^nb?wsPH ~ IO^Mq (Bate & Burkert 1997) where 
M Jeans is a Jcans mass, A^nb = 50 is the number of neighbor par- 
ticles used in the SPH calculation and mspe is the mass of an 
SPH particle. Our effective mass resolution Mjeans ~ lO^M© 
corresponds to the mass of baryonic substructure such as gi- 
ant molecular clouds and GCs in galaxies. In order to perform 
such high resolution simulation, we focus on the formation of 
a small galaxy (with total mass of lO'^M©). 

This paper is organized as follows: We describe numerical 
method and models in ^ In O we demonstrate the effect 
of low-temperature cooling (T < 10^ K) for the formation of 
clumps. In ^ we discuss the tidal evolution of clumps in the 
galactic halo. Summary and discussion are given in gS] 

2. MODELS AND NUMERICAL METHODS 

Our numerical tech nique is based on a standard A^- 
body/SPH method (e.g..' Hernquist & Katzll98 9'). where self- 
gravity, radiative cooling and star formation from dense gas 
clouds are taken into account with a high spatial resolution 
(~ 100 pc). In SPH simulations, the finer mass resolution en- 
ables us to resolve denser and colder structures. In fact, our 
particle mass (Msph ~ IO^Mq) allows us to follow the cold 
gas phase below T = 10"^ K. We thus consider gas cooling in 
the range between 10 K and 10** K. We review our models 
and numerical scheme in this section. 

2.1. Initial Conditions and Resolutions 

We study galaxy formation from a top-hat density peak. 
Additional small-scale density perturbations are added using 
GRAFIC2 in the COSMICS package ( Bertschinger 200 1). 
We use the following set of cosmological parameters: i7o = 1, 
nx = 0.0, = 0.1, h = 0.5 (the unit is lOOkm/s/Mpc), and 
(Tg = 0.63. The total mass of the sphere is ^ lO'^M©. 

The initial redshift of the simulation, Zini, is ^ 87. We im- 
pose an overdensity which corresponds to the collapse epoch 
Zc ^ 3, and the turnaround epoch Zt ~ 5.5. Since the limited 
volume cannot represent a realistic tidal field, matters in the 
region cannot get angular momentum from the external tidal 
fields. We thus initially add a rigid rotation for matters in the 
simulation sphere. The spin parameter. A, is 0.05, which is a 
tvpica l value derived fro m cosmological A^-body simulations 
(e.g.. Warr en et aljT99l . while we cannot follow the angular 
momentum evolution correctly (White 1984). 

The total number of the particles is A^dm = A^sph = 1005600. 
The corresponding masses of particles are niuu = 1.0 x 
IO'^Mq for dark matter and mspH = 1.12 x IO^Mq for bary- 



TABLE 1 
Cosmological parameters . 



rim Sl;^ h Q-g 

LP 0.0 0.1 0.5 0.63 

TABLE 2 
Parameters of the Halo. 



LOx 10'°Mq 87 5.5 3.0 0.05 
Note. Hq = lOO/i km/s/Mpc, where Hq is a current Hubble constant, erg is the 
rms density fluctuations averaged in spheres of 8 /r' Mpc radius at z = 0. 



TABLE 3 
Numerical Resolution 





particle number 


particle mass 


gravitational softening length 


DM 


1005600 


1.0 X 10'*Mq 


108 pc 


SPH 


1005600 


1.12 X IO^Mq 


52 pc 



onic particles. The gravitational softening lengths are comov- 
ingly evolved from the beginning of the simulation to z = 10 
(Governato et al. 2004), and they are fixed at 52 pc for bary- 
onic particles and 108 pc for DM particles from z = 10 to z = 0. 

Although the initial conditions we use are artificial and 
the adopted background cosmology is not a recently favored 
A-dominated universe, the background cosmology is not re- 
ally important for the semi-cosmological simulations as far as 
the magnitudes of the additional small-scale perturbations are 
reasonable. Moreover at high-redshift (z > 2) where we study 
the evolution of small galaxies in this paper, the Einstein-de 
Sitter universe is still a good approximation. 

The parameter set of this simulation is summarized in Ta- 
ble ^ (the cosmological parameters). Table |2] (the parameters 
of the halo), and Tabled (the number of particles, mass reso- 
lutions and gravitational softening lengths of particles). 

2.2. Tree+GRAPE N-body/SPH code 

We accelerate gravity calculations by combining two fast 
gravitational solvers, i.e., a special purpose computer GRAv- 
ity PipE (GRAPE; Sugimoto et al. 1990) and Tree algorithm 
(e.g., Appel 1985: Barnes & Hut.19 86). The combination of 
GRAPE and Tree (hereafter TreeH-GRAPE) is firstly proposed 
by Makino ( 1991 ) and used in astronomical simulations (e.g., 
Fukushige & Makino 200 i|). 

Parameters used in our TreeH-GRAPE algorithm are as fol- 
lows: the critical number of particles to share the same inter- 
action list, He, is 5000 (see Kawai et al. 2000). The opening 
angle for Tree, 9, is set to 0.5 and we only use the monopole 
moment. The calculation speed of gravitational forces using 
our TreeH-GRAPE code is about dozens times faster than when 
we solely use Tree or GRAPE in tests on a system with a 
GRAPE and a host. 

Gas dynamics is solved by the SPH method (e.g., iLucvl 
1977; Gingold & Monaghan 1977), in which shocks are han- 
dled via a standard artificial viscos ity. Viscosity parameters 
area= 1.0,/3 = 2.0and7/ = 0.1 ( Monaghan & Gingold 1983.) . 
The number of neighboring particles for each SPH particle is 
set to A^nb = 50. We adopt a shear-reduced technique in the 
artificial viscosity (Balsara 1995) . We solve an energy equa- 
tion of the ideal gas, 7 = 5/3, with the radiative cooling and 
the inverse Compton cooling. We assume that the gas has a 
primordial abundance ofX = 0.76 and Y = 0.24, where X and 
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Y are the mass fractions of hydrogen and helium. The mean 
molecular weight is fixed at /i = 0.59. 

We apply the Tree+GRAPE method for neighbor search 
to reduce the computational cost. We first generate a poten- 
tial neighbor list using Tree, and search neighboring particles 
from the list with GRAPE. We adopt the reordering method 
to reduce t he cost of data trans fer between GRAPE and host 
computer tSaitoh & Kod al20'03l) . The performance is an order 
o f magnitud e higher than Tree or GRAPE calculations. 

iBate & B urkert (1997) and "BatejtaD 12003) pointed out 
that, in the SPH simulations, the Jeans instability cannot 
be calculated correctly for masses less than 1.5-2A^nb'«sPH, 
where mspH and A'nb are the mass of an SPH particle and the 
number of neighbor particles, respectively. In our simula- 
tions, the thermal temperatures of gas particles are restricted 
to retain a Jeans mass larger than 2A'nb'WsPH, that is 



Mjean,s--/OX Aj, 

_^-3/2^-I/2^3>2A?„b'«SPH, 



(1) 
(2) 



where Ay, G, c,, and p are the Jeans length, gravitational 
constant, sound speed, and local density of gas, respec- 
tively. Equivalently, the gas temperature should be greater 



than Tjninip), 



Tmin(p) ■■ 



2/imp Gpi/3(2A?„bmsPH)'/^ 



3A:r 



7(7-1) 



(3) 



where nip is the mass of a proton and is the Boltzmann 
constant. In Figure^ we show critical lines of Jeans limits for 
several masses of SPH particles, which is defined by Eq.Q- 
In our model, mspH is 1.12 x 10^ Mq, therefore the minimum 
temperature of the gas with number density riH^ 100 cm"-' is 
- 500 K. 

This high resolution in our simulations allows us to cal- 
culate the cold gas phase {T <C 10"* K), which is a poten- 
tial site of star formation. Previous numerical simulations 
of galaxy formation did not have sufficient resolution to cal- 
culate such a cold phase. We adopt the cooling curve of 
Epaans & Norman ( 1997), which is used in 2-D and 3-D sim- 
ulations of the ISM on a galactic scale (e.g., Wada & Norman 
UoOl; Wada 2001). We also include the inverse Compton 
cooling Jlkeuchi & Ostrik er 1986). The allowable tempera- 
ture for each gas particle with density p is rmin(p) <T < 10^ 
K in this paper The radiative cooling below 10"* K plays an 
important role in the formation of clumps (see 

2.3. Star Formation 

Our star formation algorithm is similar to the one bv lKatzl 
([l992). We adopt a conversion type star formation. A gas 
particle which satisfies all the following conditions is eligible 
to be converted to a star particle: The gas particle (1) has 
high density («h > 0.1 cm"^), (2) lies in an overdense region 
(Pgas > 200 X Pbg(z)), where Pbc(z) is the cosmic background 
density at z, (3) has low temperature (T< 30;000K), and (4) 
is in a collapsing region (V • v < 0). 

The star formation rate of a gas particle is given by the 
Schmidt law, 

^=c.^, (4) 

at ^dyn 

where p* is the density of newly-bor n stars, the local star for- 
mation efficiency, c,, is 1/30 (e.g., 'Abadi et al."2003a^ and 
the dynamical time is fdyn = 1 / v^47rGp^. Then, the proba- 
bility of turning the gas particle into a star particle during a 




Fig. 1. — Phase diagram and limits of resolution. The Jeans instability 
of the gas above the lines can be calculated properly in SPH simulations. 
The solid, dashed, dashed-dotted, and dotted lines indicate Jeans mass limits, 
Mjeans = 2A'nb X ^SPH, with the mass resolutions of a single SPH particle be- 
ing MspH = W^Mq, Msph = iCMq, A/sph = 10^ Mq and Msph = 10* Mq, 
respectively. In the figure, we assume A'nb = 50, /i = 0.59, and 7 = 5/3. The 
contour map represents the volume-weighted phase diagram from the numer- 
ical simulation of the multiphase ISM in a galaxy ( Wada & Tomisaka 2005]). 



time step St is given by 



P= 1- 



^ fdyn ^ 



exp 



(5) 



We do not include the mechanical and radiative feedback from 
star formation, SNe and succeeding metal enrichment to the 
gas for simplicity. 

3. EFFECTS OF LOW TEMPERATURE COOLING 

In order to demonstrate the importance of the low temper- 
ature cooling (T < 10'^ K), we run two simulations with and 
without the cooling below T < lO** K from the same initial 
condition. We call the former simulation 'model A' and the 
latter, 'model B'. 

The snapshots of model A and B are shown in Figure|2] The 
figure shows dark matter (left panels), gas (central panels) and 
stars (right panels) at six different epochs from z= 14 to 2. 
Clearly, the dumpiness of the gas is very different between 
the two simulations. Model A has many clumps throughout 
the formation of the galaxy, whereas the distribution of gas in 
model B is smooth, because the gas cannot cool below lO'* K. 
In other words, small scale structures cannot be evolved due 
to high sound speed. The number of clumps in model Bis 
an order of magnitude fewer than that in model A (see ^4.H . 
Tablel^shows the numbers and mass fractions of clumps. The 
total masses contained in clumps do not change significantly 
between z = 4.7 and 3.3 in both models. Since massive clumps 
quickly sink and merge into the central galaxy between z = 3.3 
and 2.0, the mass fraction decreases to z = 2. 

Model A shows a variety of fine structures at z = 2, such as 
spirals, clumps, tidal tails, while model B does not show such 
features. These structures consist largely of the low tempera- 
ture gas. Gas cooling below lO'* K is important for investigat- 
ing the evolution of clumps, and consequently the evolution of 
host galaxies. In the following sections, we discuss evolution 
of the clumps based on the results of model A. 

4. TIDAL STRIPPING OF DARK HALOS AROUND CLUMPS 
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HadEl A Hadcl B 
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Fig. 2. — Effects of low temperature cooling (T < IC* K) on the clump formation in a galaxy. In model A, the gas can cool below T < Iff* K, whereas 
the cooling above T = lO** K is assumed in model B. Left, middle, and right panels in each model show the distribution of dark matter, gas, and star particles, 
respectively. Blue and red in each right panel indicate young (age < 0. IGyr) and old star particles, respectively. Brightness of each particle represents the local 
density of each pailicle. 



TABLE 4 

Number and Mass Fraction of Clumps. 



Model A Model B 

Number Fraction Number Fraction 



z = 4J 


54 


5.5% 


14 


2.5% 


z = 3.3 


87 


6.6% 


22 


4.2% 


z = 2.0 


62 


2.6% 


7 


0.12% 



'Number' and 'Fraction' indicate the number of clumps within the virial ra- 
dius of the galaxy and mass ratio of total mass of clumps to the mass of the 
galaxy, respectively. 



In this section, we investigate the evolution of clumps in 
the galactic tidal field of model A. We firstly define clumps 
and construct the evolutional history of each clump, we then 



quantitatively investigate the effect of tidal stripping for the 
dark halos around the clumps. We compare stripping effects 
on clumps in three different epochs, i.e., before the collapse 
epoch (z = 4.7), during the collapse epoch (z = 3.3) and after 
the collapse epoch (z = 2.0) to see the evolution. 

4.1. Identification of Clumps 

In order to identify clumps, we use SKID i Governato et alJ 
1997) as a clump finder, which is based on a local density 
maxima search. This algorithm groups particles by moving 
them along the density gradient to the local density peaks. 
The density field and density gradient are defined everywhere 
by smoothing each particle using the SPH-like method with 
64 neighboring particles. At each density peak, the localized 
particles are grouped into an initial particles list of the clump. 
Then, gravitationally-unbound particles are removed from the 
initial particles list of the clump. If the number of the bound 
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particles is more than a threshold number, nth = 100, we iden- 
tify this group of particles as a 'clump' . We define the posi- 
tions of each clump as its center of mass and the size of each 
clump, Tsub, as the distance of the particle furthest from the 
center 

In Figure [S] we plot circular velocity profiles of randomly 
selected four clumps. The circular velocities are defined as 

GM(< r)/r, where r is the distance from the center of the 
clumps and M(< r) is the mass within r, respectively. The 
upward-pointing arrow in each panel indicates the clump size 
which is derived from SKID, namely rjub- The rsubS are com- 
parable with the minimum points of the circular velocity pro- 
files (the open circles), i.e., rsub — K^c = min). Thereby, 
our clump finding method successfully finds the size of the 
clumps. 
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Fig. 3. — Circular velocity profiles of randomly selected four clumps at 
z = 2.0. The circular velocities are estimated by \/ GM{< r)lr, where r is 
the distance from the center of the clumps and M(< r) is the mass within r, 
respectively. We plot three times of r^^\y for r. The upward-pointing an'ows 
indicate the positions of r^uh^ and the open circles represent the the minimum 
points of the circular profiles, respectively. 



4.2. Distribution of Dark Matter and Baryons in Clumps 

Figure 13 shows the particle distributions at z = 2. The par- 
ticles in the clumps are plotted in the middle row. The three 
panels in the bottom row are residual particles from the orig- 
inal distribution shown in the top three panels. The positions 
of the three components, dark matter, gas, and starts, coincide 
with each other. However, the size of these components dif- 
fers, in a sense that gas and stars are centrally concentrated 
and the dark matter halos are extended. This is a natural con- 
sequence of dissipational nature of the gas. It is remarkable 
that some clumps do not have dark matter halos at z = 2 (see 
the first quadrant in the middle row of Figure 0}. 

4.3. Evolution of Clumps 

In order to trace the evolution of clumps, we construct evo- 
lution histories of clumps. Time evolution of each clump is 
tracked in the following way: At first, we prepare two clump 
catalogs (catalog / and / + 1) at two adjacent time snapshots. 
If a clump in / has more than a half of particles in a clump in 
;+ 1, it is regarded as the progenitor of the clump in catalog 
i+ 1. We trace clumps in our result back to z = 15, using 235 
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Fig. 4. — Demonstration of our clump finding algorithm. Snapshots of 
dark matter (left column), gas (center column) and star (right column) par- 
ticles at z = 2 are shown. From top to bottom, we show particles with three 
different criteria: all particles of each component (top row), the particles in 
the clumps which are detected by SKID (middle row) and the residual parti- 
cles (bottom row). To highlight structures in the halo, we plot the particles 
with the weighted brightness proportional to the local densities of the parti- 
cles. 



snapshots (time interval of 10^ yr), and we obtain the evolu- 
tional tracks of all clumps found at z = 2.0. 

Figure|5]shows the trajectories of a clump onx-y and x-z 
planes. This is a typical case of a clump being captured by 
the host halo. Figure |6l shows the distance from the center 
of mass of the clump to the galactic center (solid line) and 
the virial radius of the host halo /?vir (dashed line) against time 
and redshift. The clump gets into the virial radius of the host 
galaxy at z ~ 4.5. The pericenter-to-apocenter ratio of the el- 
liptical orbit is ^ 0.2, which is consistent to the mean value of 
observed galactic globular clusters (e.g., Dinescu et al. 199^ 
;the mean value is ^0.3) and the values in clusters in CDM 
simulations (e.g .. lGhigna et al.l l998: Okamoto & Habe 199^ 
iGiU et alJ200l . In our analysis shown below, we refer to this 
clump as the reference case. ^ 

4.4. Size of Clumps 

Since clumps in the galactic halo are subjected to a tidal 
force of the host galaxy, their sizes, r^ab, should correlate with 
a tidal radius which is determined by the tidal force. Assum- 
ing isotherm density profiles [p(r) oc r"^] for the host halo and 
clumps, the tidal radius for a clump at is defined by 



?"tidal(^h) = 



'«tot 
'•sub Mtot(^h) ' 



(6) 



where mtot and Mioi{r < Rh) are the total mass of the clump 
and the total mass of the host galaxy w ithin Rh, respectively 
jCihi gna et alJll99a lOkamoto & Habelfr999,) . This equation 
is derived from the balance between the tidal force due to the 
host halo and the self-gravity of the clump. 

Evolution of r^uh, '"tidal of the reference clump, and R^ is 
shown in Figure0 The clump is captured by the host galactic 

^ The clump labeled #34 in Figurel3l(bottom-left) is the reference clump. 
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Fig. 5. — Projected orbits of the reference clump. Solid and dashed Hnes 
indicate the trajectories on x-y and x-z planes, respectively. Origin of co- 
ordinates accord with the center of the galactic halo. Two circles show the 
starting points of each trajectory. 




Fig. 7. — Time evolution of the size of the reference clump (rsuh) and 
tidal radius (ru^i). Thick-solid and dashed Unes indicate rs„i, and rtyai, re- 
spectively. Distance of the clump from the galactic center is shown by the 
thin-solid line (see Fig.|6|. The dotted line indicates Rh before capture of the 
host galactic halo. 
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Fig. 6. — Distance of the reference clump from the galactic center, Rhi and 
the virial radius of the host galaxy, Ryi, , against time or redshift. Solid and 
dashed lines indicate and Ryin respectively. 



2010 7 



Redshift 
5 4 3 




1 ^ 2 3 

Time/Gyr 

Fig. 8. — Time evolution of the half mass radii of the total, baryon, 
and dark matter components of the reference clump. Thin-solid, dotted and 
dashed lines indicate r^f toi> ''hlf.b. ™d /"hif.dm. respectively. Thick-solid line 
indicates rs„5, which is already shown in FigurelTf 



halo at z ^4.5. The size of the clump is independent of the 
tidal radius before z ^ 4.5. However, it follows rtidai after 
z ^ 4.5. This suggests that once a clump is captured by the 
host halo, its radius instantaneously shrinks to the tidal radius. 

In Figure|8]we display the half mass radii for total (rhif.tot), 
baryon (rhif b), and dark matter (rhif dm) components as repre- 
sentative sizes of the distributions of matters within the clump. 
From the plot, we find that (1) rhif.b is always smaller than 
'"hif.dm , and (2) rhif,tot comes close to rhif,b with time while 
rhif,tot is almost similar to rhif.dm before the first pericenter pas- 
sage (z > 3.7). Due to the radiative cooling, the gas shrinks 
towards the center of the clump immediately and forms a core 
in there. This is the reason of the first point. The second 
point, namely rhif,tot shifts from ruf.dm to rhif,b, is explained as 
the effect of the tidal stripping. 

In Figure |9l we plot rtidai of clumps in the galactic halo 
against r^^b at z = 4.7,3.3 and 2.0. If the radii of clumps 
are determined by the tidal field, we should find a correlation 
(solid Une) in Figure|9] These figures clearly indicate that the 



sizes of clumps are determined by the local tidal field instan- 
taneously, as we have already shown in Figure in spite of 
the tidal stripping being more efficient at the pericenter This 
is because the matters off from clumps spill their local tidal 
radii after passing away from the pericenter. Indeed, some of 
the clumps show clear evidence of elongation as the tidal tails 
(see the distribution of old stars in model A in figure|2} 

4.5. Tidal Truncation of Dark Halos and Formation of Pure 
Baryonic Clumps 

Since the baryons have more concentrated distribution than 
the dark matter within a clump, tidal stripping should affect 
them differently. We will show this effect more quantitatively 
below. 

Figure[lOldisplays the mass evolution of three components 
(gas, star and the dark matter) of the reference clump. The 
mass of dark matter decreases with time. Mass loss oc- 
curs mainly at each pericenter passage once a clump is cap- 
tured. The rapid mass loss, for example at around f ^ 1 .3 Gyr 
iz ^ 3.7), corresponds to the time when the clump is pass- 




Fig. 9. — rtidal against r^ub- The black line indicates that both quantities coincide. 



ing the pericenter. The gas/stellar mass decreases/increases 
steadily owing to the star formation. However, the total mass 
of baryons (gas + star) is almost constant. This is because that 
the baryonic component tends to reside inside the tidal radius 
of the clump and is therefore less vulnerable to tidal stripping. 
Eventually, the total mass of the clump becomes one-tenth of 
its original mass. As shown in FigurefTO] the baryon mass to 
total mass ratio in the reference clump increases with time, up 
to 0.5. 

The baryon mass fractions in clumps in the galactic halo 
are plotted against in Figure^2(^) for ^ snapshot at z = 2. 
Many clumps show signs of tidal stripping (mb/wtot > 0.1). 
Since tidal force from the host galaxy is effective for clumps 
with smaller R^, mb/mtot should be larger for smaller R^. 
However, there is no clear trend between R^ and nib /^tot in 
Figure^2(a)- Suppose that the tidal stripping is the most ef- 
fectively occurs following each pericenter passage, we expect 
that mb/mtot is proportional to the inverse of Rpen. In fact. 
Figure [^(b) shows this trend. ^ Note also that no clump is 
found within Rp^^i < 1 kpc, since the dynamical friction time 
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Fig. 10. — The evolution of masses of the reference clump and nib/nitot 
of the clump. Thin-solid, dashed, dot-dashed, and dotted lines indicate the 
masses of baryon, gas, stars and dark matter in the clump, respectively. 
Thick-solid line shows the evolution of n% / nitot of the clump. 



' We adopt the current distance from the center of the halo, Rh. instead of 
ri, if a clump is in the first falling phase. 



scales of clumps within 1 kpc (< 10** yr) are smaller than the 
age of the system (^10^ yr) and clumps with small pericen- 
ters (< 1 kpc) fall into the galactic center and are tidally de- 
stroyed/disappeared. 

Another possible effect to determine mb /mtot is the eccen- 
tricity of each clump's orbit. Even if Rp^-i is the same, clumps 
with high eccentricities spend most of their time at apocenter 
(^apo), as a result they are less severely stripped than those 
with smaller eccentricities. In Figure[21 we plot mb/mtot as 
a function of Rpeii/Rapo at z = 2. As expected, there is a weak 
positive correlation between mb/mtot and Rperi/Rapo- In addi- 
tion to that, the plot also suggests that the baryon fraction is 
determined by (1) how many times the clumps pass the peri- 
center, represented by the numbers, and (2) Rpen, represented 
by the color sequences. 

Figure \^ shows distribution of mb /mtot in clumps at z = 
4.7,3.3 and 2.0. The vertical axis shows the number of 
clumps in each bin (0.1) normalized by the total number 
of clumps at each epoch. The total number of clumps is 
shown in Table |3 The fractions of baryon dominated sys- 
tems (mb/mtot > 0.5) increases with time as 30 %(z = 4.7), 32 
%(z = 3.3), and 60 %(z = 2.0). This is a natural consequence 
of the effective stripping of the dark matter halos from clumps 
by the tidal interaction with the host galaxy (see FigurefTol. 

5. SUMMARY AND DISCUSSION 
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Fig. 1 1 . — Distribution of m(,/miot in the galactic halo atz = 2.0. The left 
panel shows mj, / mtot against Rj, . On the other hand, the right panel shows 
m(,/mtot against the pericentric distance, Rp^i. Crosses in the panels repre- 
sent each clump. The open circle in the right panel is the clump which is 
preprocessed by the other halos. 
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Fig. 12. — Distribution of nib/mtot against Rpai/Rapo in tlie galactic halo 
at z = 2.0. Red, blue, orange, and purple numbers indicate different ranges 
of the pericentric distance from 1 < Rp^i/kpc < 2 to 4 < Speij/kpc < 5 and 
the numbers indicate how many times each clump passes the pericenter. We 
obtain pericenters and apocenters from the most recent period of clumps. 
Note that we does not consider clumps in this plot, which have never passed 
pericenters. and larger pericentric distances, Rpai > 5 kpc. 




Fig. 13. — Fraction of clumps with nib/nitot at three redshifts, i.e., z - 
4.7(solid line), z = 3.3(dashed line), and z = 2.0(dotted line). 



• Gravitationally-bound, small clumps with core-halo 
structure are first formed from CDM perturbations. Af- 
ter captured by the galactic halo, these clumps lose their 
dark matter halo due to the tidal field of the host galaxy, 
and result in baryon-dominated objects, such as GCs. 
By z = 2.0, almost half of the clumps (57 %) becomes 
baryon-dominated systems (mb/mtot > 0.5). 

• The simulation without radiative cooling below T = 10^ 
K (the model B) shows an order of magnitude fewer 
clumps. This result implies that treatment of the low 
temperature cooling as well as heating is extremely im- 
portant to investigate the formation of low mass objects, 
such as GCs. 



Clumps studied here have a mass range similar to that 
of globular clusters (GCs) or dwarf galaxies. Two forma- 
tion scenarios for GC formation has been proposed in lit- 
erature. One is that GCs are formed in DM minihalos 
before they infall into galactic halos (e.g., Peebles 1984 
Bromm& Clarke 2002; Wail & Pudritz 2001; Beasle v et alJ 
I2003f .Mashchenko & Sillsi.2005aib.) . The other is that GCs 
are formed from baryonic processes, such as hydrodynam- 
ical sh ocks, and thermal instability, in gala ctic halos (e.g., 
ISchwe izer 1987 ; Fall & R ees .1985; Whitmpre & SchweizeJ 
[l995nKravtsov & Gnedinj I'io'osI) . Clumps around a host 
galaxy in our simulation (see the last panel of Figure |2} orig- 
inate in hierarchically formed DM minihalos, and therefore 
our numerical results support the former scenario. 

The absence of DM halos around observed GCs JPrvod 
[^89; Moore 1996) could be the result of tidal stripping by in- 
teractions with the host galaxy. 'M ashchenko & Sillsl (l2005bl) 
studied the tidal evolution of primordial GCs with dark mat- 
ter halos in a static potential of a host galaxy. In their work, 
dark matter-less GCs form from proto-GCs with dark matter 
halos via tidal stripping. Our results also suggest that GCs 
can be formed naturally in the context of the cold dark matter 
universe. 

In taking some heating mechanisms such as SN ex- 
plosions and/or ultraviolet background heating and/or 
rei onization (e.g., LD ekel & Silk 1986; Efstathou 199: 
_mmer-Larsen et al. 199^ .Thacker & Couchmaa .20011 
usa & Umemura 2004a ij into consideration, the suppres- 
sion and/or destruction of small mass clumps, such as GCs 
and dwarfs, should occur. We plan a new series of simula- 
tions including feedback from stars, without sacrificing the 
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spacial and mass resolutions. 



In order to study the co-evolution of galaxy and its substruc- 
ture such as GC, we run high-resolution A^-body/SPH simula- 
tions of galaxy formation. The numerical resolution of our 
simulation is high enough (~ 50 pc and IO^Mq) to resolve 
the Jeans instability of small clumps whose mass is larger than 
1O^M0. With this high resolution, radiative cooling below lO** 
K can be taken into account and this low temperature cooling 
strongly affects evolution of these systems. What we found 
on the evolution of clumps is summarized below: 
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